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ABSTRACT 

Rapidly rotating Neutron Stars (NSs) in Low Mass X-ray Binaries (LMXBs) are 
thought to be interesting sources of Gravitational Waves (GWs) for current and next 
generation ground based detectors, such as Advanced LIGO and the Einstein Tele¬ 
scope. The main reason is that many of the NS in these systems appear to be spinning 
well below their Keplerian breakup frequency, and it has been suggested that torques 
associated with GW emission may be setting the observed spin period. This assump¬ 
tion has been used extensively in the literature to assess the strength of the likely 
gravitational wave signal. There is now, however, a significant amount of theoretical 
and observation work that suggests that this may not be the case, and that GW emis¬ 
sion is unlikely to be setting the spin equilibrium period in many systems. In this 
paper we take a different starting point and predict the GW signal strength for two 
physical mechanisms that are likely to be at work in LMXBs: crustal mountains due 
to thermal asymmetries and magnetically confined mountains. We find that thermal 
crustal mountains in transient LMXBs are unlikely to lead to detectable GW emis¬ 
sion, while persistent systems are good candidates for detection by Advanced LIGO 
and by the Einstein Telescope. Detection prospects are pessimistic for the magnetic 
mountain case, unless the NS has a buried magnetic field of B « 10 12 G, well above 
the typically inferred exterior dipole fields of these objects. Nevertheless, if a system 
were to be detected by a GW observatory, cyclotron resonant scattering features in 
the X-ray emission could be used to distinguish between the two different scenarios. 
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1 INTRODUCTION 

Rapidly rotating Neutron Stars (NSs) are considered an in¬ 
teresting source of Gravitational Waves (GWs) and are one 
of the main targets for current searc hes w ith ground based 
detectors, such as Virgo and LIGO (lRileal2013l) . The char¬ 
acteristic amplitude of the GW signal scales with the square 
of the rotation frequency, thus making the more rapidly 
rotating NSs ideal candidates for detection. In particular 
some of the most promising targets are likely to be accreting 
NSs in Low Mass X-ray Binaries (LMXBs). Not only are 
these NSs rotating with millisecond periods, but the process 
of accretion from the companion star can drive the growth 
of a quadrupolar deformation. Plausible mechanisms that 
may be at work are the creation of a “mountain” (i.e. any 
kind of non-axisymmetric deformation that gives rise to 
an l = m = 2 mass quadrupole) supported by the elastic 


crust (iBildsten i 1998 : Ushomirskv. Cu tler fe Bildstenl 


2003: 

Haskell. Jones & Andersson 2006: 

Johnson-McDaniel & Owen 20131) or by a solid c 

ore of 

unsta- 

1991; 

exotic matter (Owen 19951: Haskell et al. 2007|). 
ble modes of oscillation of the star (Anderssoii 

Andersson, Kokkotas &; Stergioulas 19991) and magneticallv 

supported mountains ( 

Cutler 2002|; Haskell et al. 

2008; 

Melatos & Pavne 20051: 

Privmak. Melatos & Pavne 

2011; 


Vigelius fe Melatosl 12009s ). 


LMXBs were originally invoked as a source of GWs 
to solve an observational puzzle. In an LMXB the NS is 
spun up by matter accreted from the companion via a disc. 
This is, in fact, how old NSs are thought to be recycled 
to millisecond periods and eventually p roduce a m i llisec¬ 
ond radio pulsar after accre t ion s tops iAIpaL_eL_alJ [ToS^; 
Radhakris hnan fc Srinivasanl 1 19821 : IPapitto et al. 2013a ). 






























































2 B. Haskell et al. 


One would therefore expect the NS to be spun up to its 
centrifugal break-up frequency, which is eq uation o f st ate de¬ 
pendent, but gener ally well above 1 KHz l|Cook et al.ll 199-1 : 
lHaensel et al.lll999l h This is not, however, what is observed. 
The distribution of spins in both LMXBs and millisecond ra¬ 
dio pulsars appear s to have a statistically significant cuto ff 
at around 730 Hz llChakrabartv et al .11 20031 : |Pat r un ol 1 20i (1 ) . 

It is natural to ask what physical process removes an¬ 
gular momentum from the NS and prevents it from spin¬ 
ning up further. The first and most obvious candidate is 
the interaction between the stellar magnetic field and the 
accretion disc. This p ossibility was examined in detail by 
IWhite fe Zhand (j 19971 ') who found that, at least for the data 
available at the time, this scenario would involve an unex¬ 
pected correlation between the accretion rate and magnetic 
field strength (which would also need to be higher than ex¬ 
pected). This led to the alternative suggestion that GWs 
may be providing the torque needed to balance the accretion 


torques, and set the spin equilibrium period of these system s 
dPapaloizou fe Pringlelll978l : I Wagonetlll98-4 iBildstenl) 1998h . 

A corollary of GW torque-balance is that the bright- 
est X-ray sour ces should also be the loudest GW emitters 
(iBildstenlll998lj . This describes the nearby LMXB Scorpius 
X-l, which has been the subject of a n umber of LIGO and 
Virgo searches dAbbott et all l2007al lbl: lAbadie et al.l l201ll : 
lAasi et aD 12014a ) that have led to a 90% confidence up¬ 
per limit for the gravitational wave strain of /i rms « 10“ 25 
around 150 Hz. With advanced detectors now coming on line 
there is a strong case to d evelop directed data analysis algo¬ 
rithms (lAasi et al .11201411) an d all-sky pipelines th at search 
for unknown binary systems dGoetz fe Rileslfeoill ). 

Although GW searches with initial LIGO are still 
not sensitive enough to probe the predictions of the GW 


torque balance scenario, the pro blem has been recently re¬ 
assess ed by several authors. IPatrun o. Has k ell fe D’An gelo 
(l201lll found that with current data the strong correla¬ 
tion between m agnet ic field and accretion rate found by 
IWhite fc Zhanel d 199711 is no longer needed and the measured 
spin period of most systems can be und erstood in terms 
of th e disc/magnetosphere interaction llAndersson et all 
l2005l l. Furthermore a detailed analysis of individual sys¬ 
tems shows that many of them do, in fact, appear 
to be close to a propeller phase in which the spin-up 
torque is much weaker than in standard accr etion models 
(lHaskell fe Patrundl201l1 : iFerrigno et al .1120131 ). Finally the 
measurements of spins and surface temperatures for most 
NSs in LMXBs are not consistent with theoretical predic¬ 
tions for GW emission due to an unstable r-mode (or at 
least not at a level that woul d allow for spin-equil ibrium 
due to torque balance) (Ho, Aridemson_& id tlaske]jjj201l]; 
Haskell. Degenaar fe Hoi 120121 : iMahmoodifar fe Strohmaveil 


20131 


GW torque balance supplies a useful upper limit to cal¬ 
ibrate searches. However if it is not the driving force behind 
pulsar spin evolution, it is natural to ask at what level the 
physical mechanisms mentioned above will give rise to GW 
emission, and whether it is likely to be detect ed. This ques¬ 
tion is crucial, given that IWatts et al.l (|200^ ~) showed that 
even at the torque balance level these systems would be chal¬ 
lenging to detect. In this paper we explore the non-torque- 
balance scenario in more detail. We focus on “mountains”, 
supported either by elasticity or magnetic stresses, and dis¬ 


cuss the level at which GW emission may be expected. We 
also take the discussion one step further and ask, given a 
GW detection, what constraints can be set on NS interior 
physics and how one could distinguish between the different 
mechanisms giving rise to the mountain using electromag¬ 
netic (e.g. X-ray) observations. 


2 THERMAL MOUNTAINS 


2.1 Crustal heating 

The outer, low density layers of a NS are thought to form a 
crystalline crust of ions arranged in a B ody Centred Cu bic 
BCC ) lattice [although recent work bv lKobvakov fe Pethickl 
2014) suggests that much more inhomogeneous configura¬ 
tions may be possible]. Above densities of 10 11 g/cm 3 , 
neutrons drip out of nuclei and form a superfluid in mature 
NSs with internal temperatures T < 10 9 K. At higher den¬ 
sities several phase transitions may occur, with nuclei no 
longer being spherical but fo rming rods and plates, the s o 
called “pasta” phases dLorenz, Ravenhall fe Pethickl 11970 ') . 
until at « 2 x 10 14 g/cm 3 there is a transition to a fluid of 
neutrons, protons and electrons which forms the core of the 
NS. 

In LMXBs accreted matter, composed of light elements, 
is buried by accretion and compressed to higher densities, 
where it undergoes a series of nuclear reactions such as elec- 
tron captures, neutron em ission and pycnonuclear reactions 
dHaensel fe Zdunfl3ll990l 'l . The observed cooling of transient 
LMXBs, as they enter quiescence, is consistent with a crust 
that has previously been heated by such reactions [see e.g. 
IWiinands. Degenaar fe Pagel ((20131) and references therein, 
although no t all details of the cooling processes ar e fully 
understood dDegenaar et al. 1 120131: ISchatz et al.ll2014l) l . 

Accretion asymmetries can produce asymmetries in 
composition and in heating, which in turn d eform the star 
and lead to a quadrupole dUshomirskv. Cutler fe Bildstenl 
i200dh Once the quadruple Q 22 is known the GW amplitude 
can be calculated as: 


16 /tt\1/2 GQ 22 9 2 


— 7(1) 


a) 


where G and c are the gravitational constant and the speed 
of light respectively, d is the distance to the source and 9 
is the angular frequency of the star. Note that we are con¬ 
sidering a quadrupolar Y 22 deformation, as this harmonic 
dominates GW emission. In this case GWs are emitted at 
twice the rotation frequency of the star. An approximate 
expression for the quadrupole due to asymmetric crustal 
heating from nu c lear reactions in the crust is given by 
dUshomirskv. Cutler fc Bildstenll2000h : 


Q 22 « 1.3 x 10 35 Rg 


\ 10 5 A" J 


Q 


30MeV 


g cm 


( 2 ) 


where is the stellar radius in units of 10 6 cm, 5T q is the 
quadrupolar component of the temperature variation due to 
nuclear reactions and Q is the reaction threshold energy. 
Higher threshold energies correspond to higher densities. 
In gene ral the reactions will heat the region by an amount 
dUshomirskv fe Rutledgell200ol) : 


ST 


(10 6 A') 


: C^ 1 p~ d 1 Q n AM 22 , 


( 3 ) 
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Source 

V 

d 

(M) 

At 

Ref. 


(Hz) 

(kpc) 

(10- 10 M o yr—!) 

(d) 


SAX J1808.4-3658 

401 

3.5 

4 

30 

Patruno et al. (2009) 

XTE J1751-305 

435 

7.5 

10 

10 

Miller et al. (’2003) 

XTE J1814-338 

314 

8 

2 

60 

this work 

IGR J00291+5934 

599 

5 

6 

14 

Falanga et al. (2005) 

HETE J1900.1-2455 

377 

5 

8 

3000 

Paoitto et al. (2013b) 

Aql X-l 

550 

5 

10 

30 

Giinebr, Giiver & Eksi ('20111 

Swift J1756.9-2508 

182.1 

8 

5 

10 

Krimm et al. (20071 

NGC 6440 X-2 

204.8 

8.5 

1 

4 

this work 

IGR J17511-3057 

244.9 

6.9 

6 

24 

Falanga et al. (2011) 

IGR J17498-2921 

400.9 

7.6 

6 

40 

Falanga et al. (2012) 

Swift J1749.4-2807 

518 

6.7 

2 

20 

Ferrieno et al. (20111 

EXO 0748-676 

552 

5.9 

3 

8760 

Degenaar et al. (2011) 

4U 1608-52 

620 

3.6 

20 

700 

Gierlinski & Done (20021 

KS 1731-260 

526 

7 

11 

4563 

Narita. Grindlav Sz Barret (2001) 

SAX J1750.8-2900 

601 

6.8 

4 

100 

this work 

4U 1636-536 

581 

5 

30 

pers. 

this work 

4U 1728-34 

363 

5 

5 

pers. 

Jlgrormi^al. (20111 

4U 1702-429 

329 

5.5 

23 

pers. 

this work 

4U 0614+091 

415 

3.2 

6 

pers. 

Piraino et al. (1999) 


Table 1. LMXBs for which we have obtained an estimate of the outburst duration At and average accretion rate 
(M). Where the reference column ind icates ‘this w ork’, we have used a fiducial power law index of T = 2 and 
the galactic absorption column from [Kalberla et al.l (120051) . We also list the distance d of the system and the spin 
frequency v. Sources in the top half of the table are AMXPs, while those in the bottom half are nuclear powered 
pulsars and their frequency is inferred from the frequency of burst oscillations, as explained in the text. We do not 
attempt to explicitly estimate the errors associated with these measurements. The most uncertain quantity is the 
distance, but our main conclusions on the detectability of the GW signals are unlikely to change unless there is a 
substantial error in the values above. 


where Ck is the heat capacity per baryon in units of the 
Boltzman constant ks, Pd is the pressure, in units of 10 30 
erg cm -3 , at which the reaction occurs, Q n is the heat 
per unit baryon (in MeV) deposited by the reactions and 
AM 22 is the deposited mass in units of 10 22 g. Note that ST 
in m is the total increase in temperature; only a fraction 
5T q /ST 1 is likely to be asymmetric in general and specif¬ 
ically quadrupolar. lUshomirskv. Cutler fc BildsterJ (|2000l l 
estimate that 5T q /ST ^ 0.1, but in reality the ratio is poorly 
known. 

After an accretion outburst, as the system returns to 
quiescence, the deformations are erased on the cr ust’s ther¬ 
mal timescale dBrown, Bildsten fc Rut led gel 1 199^ 1: 

Tth ~ 0.2 p 3/4 yr (4) 

If the system is in quiescence for longer than the thermal 
timescale in 0 , Q 22 is likely to be washed out and a new 
mountain is rebuilt during the next outburst. A shorter re¬ 
currence time, on the other hand, could lead to an incremen¬ 
tal accumulation of material. However, compositional asym¬ 
metries may be frozen into the crust, and not be erased 
on a thermal times cale, allowing for the mountain to be 
built incrementally dUshomirskv, Cutler fc BildsterJ l2000lh 
This scenario predicts the formation of large quadrupoles in 
all transient systems (10 38 < Q < 10 4 ° g cm 2 ), as we discuss 
in section I2~H The implied spin-down rate, in the case of four 
transient systems (SAX J1808.4-3658, XTE J1751-305, IGR 
J00291+5934 and SWIFT J1756.9-2508) is already excluded 
by measurements of the spin-down rate between outbursts 
dPatruno fc Wattsll2012h . We do not consider this scenario 
further, but note that if it were to occur in any transient 


system, the GW strain would be comparable to that of a 
persistent system. 


2.2 Maximum quadrupole 

Large stresses can break the crust, so one should 
also ask how large a mountain the star can sustain. 
This problem has b een studied by different authors in 
Newtonian physics (lUs homirs kv. Cutler &; BildsterJ l2000l : 
lHaskell. Jones & ; Anderssonl 120 061) and, more recen tly, in 
general relativity djohnson-McDaniel fc Owenl |2013| ). The 
results depend critically on the breaking strain a c of the 
crust, i.e. the average strain a = T/u that can be built up 
before the crust cracks, where p is the shear modulus of the 
crust and T the average stress. The breaking strain of a NS 
crust is not well known, but is known to be a c ~ 10 -2 for 
perfect crystals in a laboratory setting, and recent molec¬ 
ular dynamics simulati ons have shown that it may reach 
a c « 1CT 1 in NS crusts (jHorowitz fc Kadaul2009l ). The max¬ 
imum quadrupoles that can be sustained are thus of the the 
order of Q 22 ~ 10 38 — 10 39 g cm 2 for more massive stars 
(M « 2M 0 ) and Q 22 ~ 10 39 — 10 4 ° g cm 2 for less massive 
stars (M ~ 1.2Mq), depending on the exact equation of 
state. 


2.3 Gravitational radiation 

It is natural to ask, for the currently known LMXBs, how 
large a thermal mountain can grow and if it is detectable 













































4 B. Haskell et al. 


by current and next generation interferometers, such as Ad¬ 
vanced LIGO or the Einstein Telescope (ET). To answer 
these questions let us examine the LMXBs whose spins are 
known. They can be divided in two classes: the Accreting 
Millisecond X-ray Pulsars (AMXPs), which are detected as 
pulsars and can thus be timed, and the Nuclear Powered 
(NP) pulsars which do not pulsate, but exhibit quasi peri¬ 
odic oscillations in the tails of type II nuclear bursts. The fre¬ 
quency of these oscillations is a measure of the spin period, 
as confirmed by observations of burst oscill ations i n sources 
that are also detected as X-ray pulsars dPatruno fc Wattsl 
l2012h . The details of the LMXBs we use are presented in ta¬ 
bled] The other quantities listed in table [Tjare the distance 
to the source and additionally the average duration At and 
average mass accretion rate < M > during outbursts. The 
amount of mass that is accreted during an outburst can 
then be obtained as AM =< M > At, and inserted into 
equations © to calculate the temperature increase due to 
nuclear reactions in the crust. 

To calculate the average mass accretion rate for all 
sources we followed two different approaches. For AMXPs 
we used the data collected by the Rossi X-ray Timing Ex¬ 
plorer (RXTE) and recorded with t he Proportional Counter 
Array (PCA; see lJahoda et al.ll2006h . We used the Standard- 
2 data mode and extracted the 2-16 keV X-ray fluxes for all 
outbursts, following the procedure by van Straaten et al. 
(2003). The fluxes for each outburst were averaged assum¬ 
ing a fiducial spectral index (i.e. assuming that the spectral 
index remains constant during the course of the outburst 
and between different outbursts) taken from the literature. 
We then extrapolated to the 0.1-100 keV (bolometric) flux. 
We used the unabsorbed luminosity (where we take the hy¬ 
drogen absorption column Nh reported in the literature, 
see table [T]). The bolometric luminosity was then calculated 
from the distance in table |T| and the mass accretion rate 
as given by L acc = GM MR = Mr/c 2 . Here, L acc is the 
bolometric accretion luminosity, and we assumed a mass of 
M = 1.4 Mq, and a radius R = 10 km and tj is the con¬ 
version efficiency for the rest-mass into energy. After calcu¬ 
lating the average mass accretion rate for each outburst we 
selected (and reported in table]]]) the highest value obtained 
(i.e., we consider the biggest possible mountain). 

For the nuclear-powered accreting pulsars we used in¬ 
stead data from the All-Sky Monitor (ASM) onboard RXTE, 
which operated in the 1.3-12.1 keV band. We used the ASM 
rather than the PCA because all the 8 sources selected are 
either persistent sources or have long outbursts. The ASM, 
being a monitoring instrument, has a much better data cov¬ 
erage (although with lower sensitivity and a narrower energy 
band). In this case we selected the absorption column Nh 
and spectral index T from the literature (whenever avail¬ 
able) or, when no spectral analysis was available, used the 
galactic Nh and a simple power law model with spectral 
index P = 2. 

We caution that the results may suffer from systematic 
errors in both distance d and spectral index T. However 
the estimates are likely to be sufficiently accurate for our 
purposes. The main conclusions of this paper will not change 
unless there is a substantial error in our assumptions that 
can change the mass accretion rate by orders of magnitude 
(e.g., a substantial error in the distance). 


2.4 Transient sources 

In the left panel of figure Q] we show the GW strain corre¬ 
sponding to the maximum mountain that could be created 
during an outburst [equations © and ([Tjl]. assuming that 
5T q /5T = 0.1. We consider this to be a reasonable upper 
limit, as a significantly larger fraction 8T q /8T would lead to 
detectable pulsations in quiescence for some of the sources, 
as we shall see in section 1261 Note, however, that there is 
currently no physically motivated estimate for 5T q /5T, and 
the true value may be much smaller. We consider two cap¬ 
ture layers, a shallow one close to neutro n drip (where most 
of th e heat is predicted to be deposited llHaensel fe Zdunikl 
Il990l) ) with a threshold of Q = 30 MeV, and a deeper layer 
at a pressure of p = 10‘ i2 dyne/cm 2 , with a threshold en¬ 
ergy of Q = 90 MeV. All values lie below the maximum 
quadrupole that the crust can sustain. Note also that the 
increase in the quadrupole Q 22 due to the higher thresh¬ 
old energy Q is more than offset by the decrease in heating 
at higher pressures, as obtained from eq. ©. The results 
for the deep and shallow capture layers are thus very sim¬ 
ilar. The thermal timescale for the deeper capture layers 
is, however, Tth ~ 6 years. Hence a ‘deep’ mountain may 
never relax entirely in systems such as Aql X-l that have 
frequent outbursts, with recurrence times shorter than Tth- 
These systems may effectively behave as persistent sources 
for our purposes, and harbour larger mountains. 

In the left panel of figure 0 we compare our results to 
the sensitivity achieved by Advanced LIGO (assuming both 
detectors have the same sensitivity) and ET, first by assum¬ 
ing an integration time of 1 month (an average duration for 
an outburst) and then of 2 years. It is quite clear from the fig¬ 
ure that, even for a 2 year integration, most systems fall well 
below the sensitivity curve. Strain sensitivity curves for Ad¬ 
vanced LIGO and ET are respective ly taken fr om the public 
LIGO documenlQ and iHild, Chelkowski fc Freis3 (l2008li . A 
fully coherent search over time, T 0 b s , is sensitive to a strain 
of 

h « llAy/ Sn(v) /Tabs, (5) 

where S n (v) is the detector noise power spectral density, and 
the factor 11.4 accounts for a single tri al false alarm rate of 
1% and a false dismissal rate of 10% ll Abbott et al.ll2007al : 
IWatts et al.ll2008h . 

It is unlikely that transient systems will be strong 
enough sources for Advanced LIGO, but they are promising 
sources for ET. T his is essentially the same conclusion of 
IWatts et al] (l200Sl4 . who considered emission at the torque 
balance level, which is higher than the strai n we calcu¬ 
l ate [n ote that both our estimates and those of lWatts et al] 
(l2008h assume mountains that are smaller than the max¬ 
imum that the crust can sustain before breaking], A few 
systems appear to be close to the threshold for detection. 
However these systems are unlikely to be good targets for 
upcoming GW searches, as they have all just entered qui¬ 
escence after long outbursts, during which large amounts of 
mass were accreted and the crust was heated considerably. 
The mountain is currently relaxing on a timescale r t h and 
the recurrence time between accretion outbursts is likely to 


1 https://dcc.ligo.org/LIGO-T0900288/public 
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Figure 1. Gravitational wave strain vs frequency for mountains in AMXPs and NP pulsars. In the left panel we show transient sources, 
for which the mountain is the largest that can be created during an outburst, both in the case of a shallow (Q = 30 MeV) and of a deep 
(Q = 90 MeV) capture layer. In the right panel we show the persistent sources, for which we assume both the maximum mountain the 
crust can sustain (crosses), and a mountain that would give spin equilibrium from torque balance (solid triangles). The bars indicate the 
range given by uncertainties on the breaking strain, as described in the text. 


be long. It is thus probable that they will not be ‘on’ as con¬ 
tinuous GW sources during Advanced LIGO observations. 


2.5 Persistent sources 


For the persistently accreting sources the situation is dif¬ 
ferent. We assume that ongoing accretion builds the largest 
mountain that can be sustained. We take the quadrupole to 
be in the range 10 38 g cm 2 < Q 22 < 10 4O g cm 2 , to account 
for the un certainty in mass and equation of state, as es¬ 
timated by I Johnson-McDaniel fe Owenl (120131 ). The results 
are shown in the right hand panel of figure Q] The error 
bars account for the range discussed above. We also present 
the torque balance upper limits on Q 22 , as in Watts et al.1 
(120081) . The results for the maximum mountain comfortably 
exceed the torque balance limits. If accretion is ongoing, the 
quadrupole can thus become larger than the value needed 
for torque balance. In this scenario there is thus a net spin- 
down torque due to GW emission, and the prediction for the 
spin-down rate is: 


v 


6 x 1CT 13 


V 

500Hz ) 


( Q 22 \ 
V 10 38 g cm 2 J 


Hz/s, 


( 6 ) 


(where we have assumed a moment of inertia I = 10 45 g 
cm 2 for the star). Such spin-down is sufficiently strong to 
be detectable with current instrumentation. However, none 
of the persistent sources considered here have ever shown 
accretion powered pulsations that would allow us to test 
this prediction. Continued deep searches for pulsations from 
these objects is thus of significant importance for GW sci¬ 
ence. 

Another issue to consider is the amount of internal heat¬ 
ing required to sustain a large quadrupole. Rearranging © 
we see that, for a fiducial star of radius R = 12 km, one has: 


5T q 


3 x 10 7 


( Q 22 \ 

V 10 38 g cm 2 J 


f 30MeV 

V Q 


3 

K, 


(7) 


For hot sources with internal temperatures T = 10 s K, high 
values of the quadrupole (around Q 22 ~ 10 39 g cm 2 ) would 


require 5T q /T > 1, even for deeper capture layers. Such 
high values of 5T q /T > 1 would lead to pulsations in qui¬ 
escence at a level that is not observed. However for lower 
values of Q 22 , deeper capture layers and higher activation 
energy Q, the temperature perturbation is 5T q /T ^ 0.1. 
During accretion outbursts the resulting perturbations to 
the luminosity are 8Li,oi 5. 1 0 32 ergs~ 1 , and are not visible 
llUshomirskv, Cutler fe Bildstenll2000l ). as the emission is at 
much higher levels (Lboi = L acc ~ 10 35 —10 37 ergs -1 ). How¬ 
ever such levels of heating can make the quiescent flux vary, 
as we argue below. 

We can summarise the discussion above by asking what 
a GW detection implies for deep crustal heating. We use 
equations © and © to represent the sensitivity curve of 
Advanced Ligo and ET in terms of an equivalent quadrupo- 
lar temperature deformation ST q , as shown in figure [2] We 
consider two fiducial systems at a distance d = 5 kpc: 
a system that undergoes shorter outbursts and is colder 
(T = 5 x 10' K), for which we integrate the GW signal over 
the fiducial duration of the outburst (1 month); and a hotter 
(T = 5 x 10 8 K), persistent system for which we integrate 
the GW signal over a 2 year period. Figure [2] shows that Ad¬ 
vanced LIGO and ET will probe the ST q /T ~ 0.1 regime, 
with ET probing the possibly more realistic 8T q /T ^ 0.01 
regime. This is also the order of magnitude of the perturba¬ 
tions expected in the quiescent flux (see section 121)1) . which 
may be detectable with future X-ray satellites such as the 
Large Observatory for X-ray Timing (LOFT) or the Neutron 
star Interior Composition Explorer (NICER). 

I 11 the analysis above we make many approxima¬ 
tions. First and foremost we only consider two capture 
layers in the stars. In reality all layers contribute to 
Q 22 , le ading to larger q uadrupoles than those discussed 
above ((Ushomirskv. Cutler fe Bildstenl I 2 OO 0 I ). For shorter 
outbursts the reduced heating at higher densities offsets the 
higher quadrupole Q 22 in those regions. The reactions that 
deposit the most heat thus dominate, independently of den¬ 
sity. We are accounting for what is considered to be the mos t 
important layer at neutron drip dHaensel fe ZdunfkH l990). 











































5 Tq / T 



GW frequency (Hz) 


Figure 2. Sensitivity of current and next generation GW detec¬ 
tors to gravitational waves sourced by quadrupolar temperature 
deformations ST q in deep (Q = 90 MeV) and shallow (Q = 30 
MeV) layers of the NS crust. The gravitational wave strain is 
expressed in terms of the temperature perturbations 5T q /T that 
give rise to the mountain, as described in the text. Both 2nd 
and 3rd generation detectors will probe regimes of physical in¬ 
terest, with Advanced LIGO probing the 5T q /T « 0.1 regime, 
and ET the ST q /T « 0.01 regime. We show the sensitivity of 
Advanced LIGO both for a 1 month integration, Q = 30 MeV, 
and background temperatures T = 5 x 10 7 K (dotted curve), 
corresponding to the case of a short outburst, and for a 2 year 
integration, Q = 90 MeV and a background temperature of 
T = 5 x 10 8 K (dashed curve), more appropriate for a persistent 
system. Similarly we show the sensitivity of ET for T = 5 x 10 7 K, 
Q = 30 MeV, and a 1 month integration (dot-dashed curve) and 
T = 5x 10 8 K, Q = 90 MeV, and a 2 year integration (solid line 
curve). The region enclosed by the red box is that most relevant 
for LMXBs. 


where T is the temperature, C v the heat capacity, K the con¬ 
ductivity, p the density and e = e u — eh, with the neutrino 
emissivity and eh the energy deposition rate. To simplify our 
treatment and make a first assessment of detectability, we 
will use an n = 1 polytrope for the equation of state, and 
analytic expressions for t he contribution of ele ctrons in the 
crust to the con ductivity iFlowers & Itohlll98ll ) and specific 
heat (|Maxwelllll979l ), from which one obtains: 

K = 10 “(low^?) 7 (®i<) ,rg/(cmBK) ' (9) 

- 3 - 72 x 1(l " (£) (iPk) erg/<cm> K) ' (10) 

where po is the nuclear saturation density. We set 
th = 0 (note this is true for the background, but 

for the perturbations we will have 5th ^ 0), and 

for th e neutri n o emissi v ity we app rox imate the results 
of iHaensel. Kaminkerfe Yakovlev! dl996lf for vv electron 
Brehmstrahlung as: 

e„ = 6.46 x 10 18 ( —A-t) ( 7 ^ 7 ) cm 2 /s 3 (11) 

\ 10 12 g/cm 3 / V 10 9 K y ' 

At the boundary with the core we assume a constant temper¬ 
ature and for the outer boundary we assume that the emis¬ 
sion from the surface is thermal, i.e. —KX7T = (R 2 /R 2 )aT^, 
with a the Stefan-Boltzman constant. The stellar radius is 
R, and R » is the radius at which we fix the outer bound¬ 
ary of our numerical grid. The surface temperature T s at R 
is then obtain ed from the temperature T at f?„ using the 
prescription of iGudmundsson, Pethick fe Epstein! (1198311 : 




( 12 ) 


so unless there is significantly more heating in deeper lay¬ 
ers that previous calculations have not accounted for, it is 
unlikely that our results severely underestimate Qi 2 - In gen¬ 
eral the result of our analysis is that thermal mountains on 
NSs in transient LMXBs are likely to be very challenging to 
detect, even with third generation detectors. Persistent sys¬ 
tems, however, offer a promising target and electromagnetic 
observations may allow further constraints on the physics of 
the system, as we shall see in section [3711 


2.6 X-ray flux variations 

We now focus on the observable X-ray flux variations in¬ 
duced by a thermal perturbation due to a mountain in the 
crust. As already discussed we restrict our attention to per¬ 
turbations of the thermal quiescent emission. We thus as¬ 
sume that a mountain has been created during an accretion 
outburst, and study how the associated thermal perturba¬ 
tions evolve as the system returns to quiescence. To un¬ 
derstand how the surface flux is affected we consider the 
quadrupolar flux variations as linear perturbations on a 
spherically symmetric background. We start by obtaining 
the spherical background model for the temperature profile 
from the Newtonian heat transport equations in the crust: 

C v 77 - = V (KVT) — pt, ( 8 ) 


with g i 4 the gravitational acceleration in units of 10 14 cm 
s' 2 . Note that one can model th e compositio n of the outer 
laye rs in more d etail (see e.g. lHaskell. Deeenaar fe Holl2012l 
and iMahmoodifar fe Strohmaveill2013lf . However, given the 
many simplifying assumptions and the uncertainties associ¬ 
ated with the measurements in tablcQ] we use the expression 
in m, as it is unlikely to be the main source of error in our 
analysis. 

We obtain our background model by specifying a core 
temperature at the inner boundary (the crust/core interface) 
and evolving equation m until we obtain an equilibrium. 
We are now ready to evolve the quadrupolar temperature 
perturbations due to the mountain on this background. The 
evolut ion equatio n f or an l = m = 2 pertur bation takes the 
form (|hTshomirskv. Cutler fe Bildstenll2000l l : 


85T q 

Cv ~dT 


1A (r^K^h) 

r 2 dr \ dr J 


1(1 + 1 ) 


K5T q 


+ 


f 5K 
\K 



+ Fq 


dr 



(13) 


with Fq = —KS7T the background flux obtained 
from the equilibrium solution of ©■ We obtain 5K 
from equation (|9]) with the conditi on 5p = 0, as in 
llUshomirskv. Cutler fc Bildstenl l200d h At the boundary 
with the core we assume that 5T q = 0, while at the outer 
boundary we perturb the thermal flux condition, so that one 
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has SF = 4(R 2 /r 2 )aTg 5T°, with 




T 


We take <5e = 8th and assume 8th to be due to quadrupolar 
energy deposition in the capture layers. For the deep cap¬ 
ture layer, we specify an energy deposition term 8th with 
a Gaussian radial profile, located at a pressure of P = 10 32 
dyne/cm 2 , and with a half width of 5 m for the deep capture 
layer, and at P = 10 3 ° dyne/cm 2 and with a half width of 1 
m for the shallow layer. Evolving equation m we find that, 
as the problem is linear in the perturbations, to a very good 
approximation the following relations hold: 


8F_ 
Fq 
8F 
Fq 


1.29 -j?-, (deep layer) 
1.48 —(shallow layer) 


(15) 

(16) 


with very little dependence on the chosen background tem¬ 
perature T. We remind the reader that we are normalising 
to the quiescent (thermal) flux Fq obtained from the equi¬ 
librium solution of (|5|). 

In quiescence the quadrupolar temperature perturba¬ 
tions associated with a mountain and gravitational wave 
emission [equations ([2]) and ([3])] thus perturb the X-ray flux 
from the surface [equations (1151) and (1161) 1. and as the star 
rotates this leads to pulsations at twice the rotation fre¬ 
quency (i.e. the same frequency as the gravitational waves). 
In figure [3] we show the sensitivity curve for Advanced LIGO 
and ET in terms of an equivalent pulsed fraction of the X- 
ray flux. We can see that if it is possible to integrate the 
signal for 2 years (physically this corresponds to a capture 
layer deep enough that Tt.h 2 years, and the mountain 
is not dissipated significantly during the observation), both 
Advanced LIGO and ET can probe an interesting region of 
parameter space, with SF/Fq < 0.01. 


3 MAGNETIC MOUNTAINS 


3.1 Hydromagnetic evolution 


Accretion not only perturbs the structure of the star by 
affecting nuclear reactions in the crust, but it also de¬ 
forms the stellar magnetic held. As matter is accreted 
and spreads towards the equator it drags the held with 
it, and compresses it. This leads to a l ocally strong hel d 
that c an sustain a “magnetic” mountain ([Paynefc Melatod 
120041 ; iMelatos fe Pavnel 120051 : IVieelius fe Melatosl 2009bh . 


This can lead to much larger deformations than those due to 
the overall background magnetic held, even if the internal 
toroidal component of the held is much stronger than th e 
inferred external magnetic dipole (ICiolh fe Rezzollal 120131 . 
Recent calculations have shown that for realistic equations 
of s tate t h e mountain cou ld lea d to a detectable GW sig¬ 
nal (IPrivmak, Melatos fe Pavnell 20 lih . Note also that mag¬ 
netic mountains are not sustained by crustal rigidity and 
the resulting quadrupole can thus be larger than the value 
required to crack the crust. 

One of the main differences with respect to thermal 
mountains is that the time-scale on which a magnetic de¬ 
formation relaxes, after an outburst, is not the thermal 



GW frequency (Hz) 

Figure 3. The pulsed fraction 5F/Fq corresponding to a GW 
detection at threshold [as obtained from equation Q ] for Ad¬ 
vanced LIGO and ET. For deep capture layers (Q = 90 MeV) 
and a background temperature of T = 5 X 10 8 K, we show the 
results for an integration time of 2 years (dashed curve for Ad¬ 
vanced LIGO, solid curve for ET). Physically this is due to the 
fact that in the deep crust r t h > 2 years, and the mountain will 
thus not relax significantly during the observation. For shallow 
capture layers (Q = 30 MeV) and T = 5 X 10 7 K, we show the 
results for an integration time of 1 months (dotted curve for Ad¬ 
vanced LIGO, dot-dashed curve for ET). The region enclosed by 
the red box is the region of interest for LMXBs. 


timescale r t h, but the slower Ohmic diss ipation timescale 
T 0 > 10 8 years (IVigelius fe Melatosl i2009bl) . Hence a moun¬ 
tain forms gradually over several outbursts. Grad-Shafranov 
calculations indicate that the hydromagnetic structure of 
a mountain conforms to a single-parameter family of so¬ 
lutions which, once the size of the accreting polar cap is 
fixed, are function only of the mass accreted over the sys¬ 
tems lifetime, M a . This suggests that magnetic mountains 
can be treated as the persistent sources of the previous 
section. The main difference is that the quadrupole does 
not depend on crustal rigidity, but on the magnetic field 
strength when accretion begins, which we denote £>« (note 
that this is different from, and generally lower than, the 
expected NS magnetic field at bir th, as obta ined for m po p- 
ulation synthesis models rtFaucher- Gigucrc & Kaspi 2006JH, 
the initial field structu re, and M „. (IMelatos fe Pavnel [29051: 
IPavne fe Melatosl 120041 : IPrivmak. Melatos fe Pavnel 2011). 


As more mass is accreted the external dipolar com¬ 
ponent of the field, B ex t, is quenched according to 
IShibazaki et ali (1 19891) : 

-1 


B ex t — B * ( 1 + 


Ma 

M c 


and the mass quadrupole is given by 


Q 2 


10 4t> A 


M a 

M e 


1 + 


Ma 

M c 


g cm 


(17) 


(18) 


where f w 1 is a geometric factor that depends on the 
equa tion of state and accretion geometry (IMelatos fc Pavng 
12005I ). while M c is the critical amount of accreted mat¬ 
ter at which the mechanism saturates, w hic h also de- 
pends on the equation of state (IPrivmak. Melatos fe Pavnel 
















































Figure 4. GW strain vs frequency for the systems in table[l] for two magnetic mountain scenarios. In the left panel we show the strain 
that can be achieved assuming that the magnetic mountain does not decay between outbursts, for two values of the magnetic field, 
B * = 10 10 G and B t = 10 12 G. In the right panel we consider the scenario in which the mountain decays between outbursts. Detection 
will be very challenging for both Advanced LIGO and ET, unless B t m 10 12 G before the onset of accretion. 


1201 ll '). The estimates above are valid to leading order 
in M a /M c ; for M a ~ M c they are no longer accurate 
and numerical solutions are necessary. General relations 
for th e cri tical mass were d erived by iMelatos fc Pavn3 
(1200511 and IPavne fc Melatosl (120041 ') for isothermal moun¬ 
tains, while for more realistic equations of sta te [models 
C and E of IPrivmak. Melatos fc Pavnel (l 201 ll )l. one has 
M c r; 10" 7 (B,/10 12 G) 4 / 3 M q . In the regime M a > M a , 
both relations are expected to deviate significantly from the 
simple estimates above in G3 and m- Numerical simula¬ 
tions cannot probe this regime; instead one finds that, for 
M a > M c one has 0.01 < B ex t/B * ^0.1 and quadrupoles 
in the range 10 3 ' < Q 22 < 10 38 g cm 2 , for an initial field 
of B « = 10 12,5 G. Note, however, that the main difficulty 
in pushing the simulations to M a > 10 M c is numerical. 
The only firm upper limit on the suppression of the ex¬ 
ternal dipole field come from Ohmic diffusion, which lim- 
its the burial of the field at a level of B ex t/B * ~ 1CP 4 
dVigelius fc Melatosll2009bh . 


3.2 Pre-accretion magnetic field 

What limits can we set on B*, the strength of the magnetic 
field at the onset of accretion? Observational constraints can 
be obtained from measureme nts of the spin down bet ween 
outbursts for 4 systems [see IPatruno fc Wattsl J 2012 h and 
references therein], which are consistent with B ex t ss 10 s 
G. The magnetic fields inferred for millisecond radio pulsars 
are also in the range B ex t ~ 10 8 G. Furthermore observa¬ 
tions of a slow (11 Hz) pulsar in Terzan 5, IGR J17480-2446, 
indicate that this system, which is thought to have been ac¬ 
creting for a shorter period of time than most of the LMXB 
po pulation, may have a stronger field 10 9 G < B ex t < 10 10 
G (ICavecchi et al.ll201ll ~). It is thus plausible that one starts 
with B , > 10 9 G, and that the external field is reduced to 
B ex t ~ 10 s by accretion. 

If B » > 10 11 G, and polar magnetic burial is very 
short lived, we would expect B ex t ~ B * in the millisecond 
radio pulsars [unless accretion leads to significant dissipa¬ 


tion of the field jKonar fc Bhattacharval[l997l . 119991 )1. This 
would lead to larger spin-down rates than those observed. 
On the other hand, if the field remains buried and the mag- 
net ic mountain is stab le on lon g timescales [as simulations 
by |Vfgelhjg__&^Jvlelatpd (120091;) indic ate], then the results 
of IPrivmak. Melatos fc Payne d201ll4 imply a quadrupole 
Q 22 > 10 36 (B*/10 n G ) 4 / 3 g cm 2 for M a = M c . From equa¬ 
tion ©. this gives a spin down rate v > 10 14 Hz/s for a 
500 Hz pulsar, close to the maximum spin down rates mea¬ 
sured for millisecond pulsars. High initial fields of B» > 10 11 
G would thus challenge current observations. 


3.3 Gravitational radiation 

In the left panel of figure [4] we plot GW strain vs frequency 
for magnetic mountains that do not decay between out¬ 
bursts. We consider first a scenario in which B * ~ 10 10 
G and the critical mass M c has been accreted over a sys¬ 
tem’s lifetime (i.e. M a = M c ). The GW emission is pre¬ 
dictably weak. Given the uncertainties associated with mod¬ 
elling field burial, in figure [4] we also consider the case in 
which the birth field is B * ~ 10 12 G. In this case some of 
the systems could be emitting detectable gravitational ra¬ 
diation, and a detection would provide evidence for a high 
degree of field burial. The latter scenario can be excluded in 
the three syste ms [SAX J1808.4-3658, XTE J1751-305, IGR 
J00291+5934 (IPatruno fc Watt3l20ld 'l1 for which we have a 
measured spin down between outbursts. In all cases the spin 
down rate is v « —10 -15 Hz/s and it implies an upper limit 
of Q 22 ~ 10 36 g cm 2 , if we assume that GW emission is the 
dominant spin-down mechanism. GW emission at this level, 
due to a magnetic mountain, implies B, ~ 5 x 10 10 G and 
would be unlikely to be detected, as can be seen from figure 
Q] For our models with M a = M c , B » r 5 x IQ 10 G lead s 
to Bext ~ 2.5 x 10 10 G (|Privmak. Melatos fc Pavnell20li] ). 
Such a strong dipole field would, however, lead to a greater 
than observed spin down due to magnetic dipole radiation. 
In fact, if the spin down is attributed to dipole radiation, the 
implied magnetic field is B ex t ~ 10 8 G for all these systems 
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dPatruno fc Watt3l2012l ). An upper limit of \i/\ < 2 x 10 -15 
Hz/s also exists on the spin-down rate of Swift J1756.9-2508. 
In this case the limit on the dipole field from electromag¬ 
netic spin down is of Bext <5x 10 8 G but the field needed 
to explain the spin down in terms of GWs from magnetic 
mountains is B* « 10 12 G corresponding to B ext « 5 x 10 11 
G for M a = M c in our models. It is important to note 
though that while simula tions indicate that the quadrupol e 
saturates for M a > M c (IWette, Vigelius fe Melatosl l2Q10i) . 
no such effect is observed for the decay of the exter¬ 
nal field, and the limits on evolving the field further are 
mainly numerical. One cannot thus exclude high degrees 
of field burial. In fact the harmonic content of thermo¬ 
nuclear bursts suggests that in some systems burning oc¬ 
curs in patches and is confined by locally strong an d com ¬ 


pressed magnetic fields fehattacharvva fe Strohmavei] 


Misanovic ii ^allowa^_^_Coo2ei 20ld: ICavecchi et al 


Chakrabortv fe Bhattacharvvall2012l ). 


2006; 


2011 ; 


We also analyse the scenario in which the magnetic 
mountain decays on short timescales between accretion out¬ 
bursts. Time-dependent MHD simulations show that mag¬ 
netic line tying at the stellar surface stabilises the mountain 
against interchange instabilities. Current-driven Parker-type 
instabilities do occur, but they do not disrupt the moun¬ 
tain, satur ating in a state where the quadrupole is reduced 
by < 60% llVigelius fe Melatosll2009bh . Simulations confirm 
stability up to the tearing-mode timescales but they do not 
resolve slower instabilities and modes below the grid-scale. 
Different ch oices of b ou ndary conditions can also d estabilise 
the system dMukheriee, Bhattacharva fe Mignondl2013a| [bli. 
In this scenario we take M a = At(M) for each system, and 
calculate the quadrupole from equation m■ The results 
for the predicted gravitational wave strain are shown in the 
right panel of figure [4] This scenario leads to small moun¬ 
tains and weak GW emission, that would be undetectable 
for most systems, even for ET. The only systems that would 
lead to detectable GWs are the persistent ones, if B, « 10 12 
G. 


In figure[o]we show the GW strain expressed in terms of 
an e quivalent B ex t obt ained from equa tion fl8l using model 
E of lPrivmak. Melatos fe Pavnd (1201 ll'l and M a = M c , for 
which B ex t = B*/2. We can see that Advanced LIGO is ex¬ 
pected to probe high field scenarios, with 10 11 G < B, < 
10 12 G, while ET will probe a physically more realistic sec¬ 
tion of parameter space, i.e. B * < 10 11 G. 


3.4 Distinguishing magnetic from thermal 
mountains 

An interesting question is if, given a GW detection, it would 
be possible to understand whether we are observing a ther¬ 
mal or magnetic mountain. We have already discussed the 
electromagnetic counterpart of a thermal mountain in sec¬ 
tion EU and showed that a quadrupolar deformation could 
lead to flux modulations and pulsations in quiescence at 
twice the spin frequency. The results of the previous sec¬ 
tion suggest that if a magnetic mountain were to be de¬ 
tected in a hypothetical system, such a NS would have 
a strong ’birth’ (i.e. at the onset of the LMCB phase) 
magnetic field B * ~ 10 12 G, although the external dipo¬ 
lar field may be lower, due to accretion induced magnetic 
burial. In such a circumstance cyclotron resonance scat¬ 



Figure 5. The sensitivity of Advanced LIGO and ET to a mag¬ 
netic mountain. The GW strain is expressed in terms of the mag¬ 
netic field B e x t of the star, for a fidu cial sy stem at 5 kpc and 
model E of lPrivmak, Melatos fc Pavnel (|201 lh . We take M a = M c 
and, as described in the text, one has Bext = B*/2 for these mod¬ 
els. We plot both the case of a 1 month integration (dot-dashed 
curve for ET and dotted curve for Advanced LIGO) and a 2 year 
integration (solid curve for ET and dashed curve for Advanced 
LIGO). We can see Advanced LIGO will probe high field scenar¬ 
ios, with 10 11 G < B ex t < 10 12 G, while ET will be able to probe 
fields of B ex t < 10 11 G. 


tering features should appear in the X-ray emission and 
IPrivmak. Melatos fe Laskvl (120141 ) have studied the problem 
in detail for the case of an accretion buried field. 

We repeat the analysis here for a 1.4 Mq NS with an 
ac crete d outer env elo pe descr ibed by equation of state E 
of lPrivmak. Melatos fe Laskvl d2014i ). We vary B* between 
10 11 G and 10 12 G and study the emission features for M a = 
Me- In figure [6] we show an example of the kind of spectra 
that such a setup produces. The solid line represents the 
phase averaged spectrum, while the dotted lines represent 
phase resolved spectra for two extreme rotational phases, 
u> = n/2 and uj = 37r/2. We can see that in all cases the 
energy of the first line is fairly stable, but the depth can vary 
strongly with phase, as can the shape of the higher energy 
features. A strong phase dependence of the fundamental line 
for diff erent sizes o f polar m ountains h a s als o been found by 
iMukheriee, Bhattacharva fe Mignonel d2012l ). 

Let us focus on the phase averaged spectrum. Our sim¬ 
ulations show that for B* < 10 12 G no cyclotron resonance 
scattering features are present. The results for higher field 
strengths are shown in figure 0, where we plot the dif¬ 
ference in depth between the first and second line and the 
ratio between the energies at which the lines appear, versus 
the pre-accretion magnetic field B*. The effects are small, 
but may be measurable by future X-ray observatories such 
as NICER and LOFT, which will both be capable of re- 
solving modulations of less than 1% at e nergies of « 1 keV 
(iGendrau et al.l 120121 : iFeroci et al] 120121 ). Furthermore the 
cyclotron features appear to be more pronounced in the re¬ 
gion of interest, i.e. the field strengths that would lead to 
GW emission at the Advanced LIGO and ET threshold. 
This method thus has the potential to be a good diagnostic 
for distinguishing different kinds of continuous GW emis¬ 
sion. Additionally instruments such as NICER and LOFT 
will also be able to carry out phase resolved spectroscopy, 
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Figure 7. Difference in depth D in normalised flux units (left panel) and ratio between the ene rgies (right panel) of the second and 
first cyclotron line for a 1.4Mg NS obtained with model E of lPrivmak. Melatos K- Laskv { 2014 1. for varying field strengths B t . The 
different colours represent different inclinations of the observer (/.) and of the magnetic filed axis (a) with respect to the rotational axis: 
t = 0,a = 0 (blue squares), i = 0, a = 0.5ir (green diamonds), t = 0. ott , o: = 0.5 tt (red triangles), and /, = 0.257r, a = 0.257T (orange 
circles). For B r < 10 12 the features cannot be distinguished. 



Figure 6. Examp l e of a cyclotron s pectru m, obtained with the 
code of IPrivmak, Melatos fe Laskvl ll2014h for a M = 1.4Mq 
NS described by EOS E, with the following parameters: i = 
7r/4 (observer inclination relative to the rotation axis), a = 
77/4 (inclination of the magnetic axis relative to the rotation 
axis), B t = 10 12,8 G, M a = M c = 3.01426 x 10~ 7 Mq (see 
IPrivmak. Melatos fc Laskv (12014 1 for a full description of the pa¬ 
rameters). The solid line represents the phase averaged spectrum, 
while the dashed and dotted lines represent the phase resolved 
spectra for two extreme rotational phases, lo = 7r/2 (dashed line) 
and u = 37t/ 2 (dotted line). While the energy of the lines remains 
fairly constant the depth varies significantly with phase. The flux 
is normalised to give unit peak flux. 


allowing for a much more detailed characterisation of the cy¬ 
clotron resonance scattering features in these systems, which 
can vary significantly with phase, as illustrated in figure [6] 
We stress here that no cyclotron lines have been de¬ 
tected in LMXBs containing neutron stars rotating with mil¬ 
lisecond periods, and the fields of these systems are generally 
thought to be reasonably weak (B ex t ~ 10 s G). Neverthe¬ 
less if an as yet unobserved system (e.g. a system that is 


currently in quiescence) were to become visible and emit 
detectable GWs, the presence of a cyclotron line would 
point to a magnetic mountain. Its absence, on the other 
hand, combined with the estimates in section [5] would sug¬ 
gest that the quadrupole is more likely to be due to ther¬ 
mally generated crustal mountain [although mountains in 
the core of the star are also a possibility (lHaskell et al.l 
l2007i) ]. Especially for weaker fields, however, several com¬ 
binations of orientation and inclination could lead to cy- 
clot ron resonance scatte ring fea tures not being detectable 
(see IPrivmak. Melatos fc Laskvl ( 2014 ) for an in depth dis¬ 
cussion), so their absence is inconclusive. 


4 CONCLUSIONS 


In this paper we asses the likely gravitational wave sig¬ 
nal strength and detection prospects for deformations, or 
‘mountains’ on neutron stars in LMXBs. Unlike most previ¬ 
ous work on this topic we do not assume that the grav¬ 
itational wave spin-down torque has to balance the ac¬ 
cretion induced spin-up torque, as several studies have 


indicated that this is un l ikely t o be the case for many 
sy stems (lAndersson et ah 20051: Haskell fe Patrunj l201ll : 


iPatruno. Haskell fc D’Angelol 1201 ll ) . Rather, we calculate 
the gravitational wave signal strength due to the two main 
mechanisms that have been suggested for building a moun¬ 
tain: asymmetric thermal deposition in the crust (thermal 
mountains) and magnetically confined mountains (magnetic 
mountains). We calculate the gravitational wave strain for 
both mechanisms in known LMXBs for which we can mea¬ 
sure the spin frequency, average accretion rate during out¬ 
bursts and outburst duration. 

One of the main uncertainties is the timescale on which 
the mountain is stable once accretion ceases and the system 
enters quiescence. For thermal mountains it is likely that 
the quadrupole will dissipate on a thermal timescale Tth 6 
years, leading to large mountains only in persistently accret¬ 
ing systems. I11 this scenario the GW signal strength for most 
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transient systems falls below the level that would be de¬ 
tectable by Advanced LIGO or ET. In the case of persistent 
systems, however, the mountain could be even larger than 
what is required for torque balance, if the crust is as strong 
as predicted by simulations llHorowitz fe Kadaull2009i l. This 
would not only lead to detectable gravitational waves, but 
also predicts a spin-down rate of the neutron star that could 
be measurable if accretion powered pulsations were to be 
discovered from these systems. Continued deep searches for 
pulsations from luminous LMXBs, such as Sco X-l, are thus 
co mplimentary to ongoing GW searches from these systems 
dAasi et all l2014bll and could provide crucial constraints. 

For the magnetic case simulations indicate that the 
mountain could be stable on long timescales (essentially the 
Ohmic dissipation timescale r 0 ~ 10 8 years), building up 
over multiple accretion outbursts. The size of the mountain 
is strongly dependent on the strength of the magnetic field 
when accretion begins, B*. This is not well constrained, but 
the systems we consider are old systems, in which the mag¬ 
netic field is thought to have decayed and to be weak. For 
both LMXBs and millisecond radio pulsars (that are ex¬ 
pected to form mostly from LMXBs) the inferred exterior 
field strengths are B ex t ~ 10 8 G. The exterior dipolar field 
will, however, be quenched as the magnetic field is buried 
by accretion. Our simulations suggest that the exterior field 
will be reduced by approximately two orders of magnitud e 
llPavne fe Melatoa 120041 ; IPrivmak, Melatos fe Pavng|20lih , 
but the process does not appear to saturate, and the limits 
on pushing the results further are mainly numerical. Further 
field burial is thus possible. We consider two scenarios: one 
in which B* = 10 10 G and the other in which B* = 10 12 
G. For B» ~ 10 10 G the detection prospects for magnetic 
mountains are pessimistic. For a detection with Advanced 
LIGO or ET it is necessary to have an initial magnetic field 
B» « 10 12 G. Although this appears unlikely for currently 
observed LMXB systems, for w hich the evidence sug gests 
weakly magnetised neutron stars dD ’Angelo et al.|[2014l i , the 
process of magnetic burial is still not well understood, and 
such high values of the background field cannot be excluded. 

Finally, it is interesting to note that if a mountain is 
detected by LIGO or ET, it could be possible to distinguish 
between a thermal or a magnetic mountain. For the rela¬ 
tively high values of the magnetic field B* ~ 10 12 G that 
make the magnetic mountain detectable one would, in fact, 
expect phase-dependent and non-trivial cyclotron resonance 
scattering features to be present in the X-ray spectrum. We 
calculate examples of such features and show that they could 
be detected by future X-ray observatories, such as LOFT 
or NICER.. A detection of a gravitational wave signal com¬ 
bined with a detection of cyclotron features would provide 
a strong direct indication of a magnetic mountain and of a 
large buried magnetic field. 
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